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INDEPENDENT ANALYSIS OF THE ORBITS OF PIONEER 10 

AND 11 



VIKTOR T. TOTH* 

Independently developed orbit determination software is used to analyze the orbits of 
Pioneer 10 and 11 using Doppler data. The analysis takes into account the gravitational 
fields of the Sun and planets using the latest JPL ephemerides, accurate station locations, 
signal propagation delays (e.g., the Shapiro delay, atmospheric effects), the spacecrafts' 
spin, and maneuvers. New to this analysis is the ability to utilize telemetry data for spin, 
maneuvers, and other on-board systematic effects. Using data that was analyzed in prior 
JPL studie»H2l the anomalous acceleration of the two spacecraft is confirmed. We are 
also able to put limits on any secondary acceleration (i.e., jerk) terms. The tools that 
were developed will be used in the upcoming analysis of recently recovered Pioneer 10 
and 11 Doppler data files. 

Keywords: Fundamental physics; Pioneer anomaly; gravitation; solar system dynamics; 
deep space navigation; radio science; telemetry. 



1. Introduction 

As first reported in 199cm, and later confirmed in 2002^1, the Pioneer 10 and 11 
spacecraft have experienced a small, approximately constant, sunward-pointing ac- 
celeration during the later interplanetary cruise phase of their missions. This ac- 
celeration is approximately 8.74 x 10 -10 m/s 2 in magnitude; within the achievable 
margin of error, its value remains constant during the time period spanned by the 
data sets (11.5 years for Pioneer 10, 3.75 years for Pioneer 11) that were investigated. 

Although a satisfactory explanation is yet to be found, the authors of the 
two studiea^l stressed that the most likely cause of the anomalous acceleration 
is an unaccounted for on-board systematic force, such as a fuel leak, outgassing, 
or anisotropic thermal radiation. Although an attempt was made to account for 
these forces, due to the lack of available information no definitive conclusions were 
reached. 

Recently^ we successfully recovered the entire mission record of Pioneers 10 and 
11 in the form of telemetry data files containing Master Data Records (MDRs.) The 
ultimate goal of the effort described in this paper is to develop orbital analysis tools 
that can estimate the thermal emissions from the spacecraft and the associated 
recoil force using telemetry. Before that effort can proceed, however, it is necessary 
to validate our tools by re-analyzing previously studied data sets. 

We begin this paper by briefly describing the Pioneer spacecraft and missions in 
Section [2] In Section [3l we describe an independently developed approach for orbit 
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determination. Section [4] presents results: the anomalous acceleration is confirmed, 
validating the tools presented here and demonstrating that they can be used to 
achieve sufficient precision to study the Pioneer spacecraft. We also characterize the 
extent to which a jerk term may be present in the Doppler data. Lastly, Section [5] 
describes our plans to use these tools and related expertise to explore recently 
recovered Doppler and telemetry data, the goal being a satisfactory and complete 
description of the twin Pioneer missions in their entirety, from launch (or earliest 
Doppler data point available) until the end of mission. 



2. The Pioneer 10 and 11 spacecraft 




Fig. 1. The Pioneer spacecraft. 



2.1. The Pioneer spacecraft 

Launched in 1972 (Pioneer 10) and 1973 (Pioneer 11), these two Pioneer spacecraft 
were the first to traverse the asteroid belt, fly by Jupiter, Saturn (Pioneer 11 only), 
and explore the outer solar system. Originally designed to operate for at least 18 
months, the spacecraft far exceeded expectations: Pioneer 11 remained functional 
until 1995, whereas Pioneer 10 sent its last signal more than 30 years after its launch, 
in 2003. 
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These two Pioneer spacecraft remain the most precisely navigated deep space 
vehicles to date. This is due to several factors, most notable among them is the fact 
that the spacecraft are spin stabilized: therefore, they require relatively few attitude 
correction maneuvers, which means that most of the time, they fly undisturbed, 
their trajectories governed solely by the laws of celestial mechanics. 

In appearance (Figure[T]), the Pioneer spacecraft are dominated by a 2.74 m high 
gain antenna (HGA). The spacecraft body is situated behind the HGA, and consists 
of a main compartment housing most spacecraft systems, and a science compartment 
that contains science instruments. Two approximately 3 m long booms each hold a 
pair of radioisotope thermoelectric generators (RTGs) (4 in total). A third boom, 
approximately 6 meters in length, holds a magnetometer. The total (wet) mass of 
the Pioneer spacecraft is ~260 kg, of which ~30 kg is propellant and pressurant. 

The main body of the spacecraft is well insulated, designed to maintain ac- 
ceptable operating temperatures for all on-board systems both in the near-Earth 
environment, exposed to the Sun, and also far from the Sun in deep space. A pas- 
sive louver system with springs and bimetallic actuators, radiatively coupled to the 
spacecrafts' electronics platforms, ensured that excess heat was radiated away in 
space when the spacecraft was warm, but also that heat was retained when the 
spacecraft was far from the Sun. 

The Pioneer spacecraft communicated with the Deep Space Network using S- 
band microwave transmissions. The nominal uplink frequency was approximately 2.1 
GHz, the nominal downlink frequency 2.29 GHz. The spacecraft provided no range 
observable. The spacecraft also lacked a precision oscillator on board; however, their 
transmitter could operate in coherent mode. In this mode, the transmitter locked on 
to the frequency of the received signal, and the transmitted signal had a frequency 
that was an exact rational multiple (240/221) of the received signal's frequency. In 
this mode, the accuracy of the received signal was limited only by the stability of 
the ground-based transmitter, which allowed precision Doppler observations to be 
made. 

2.2. The Pioneer anomaly 

The Pioneer anomaly is best described as our apparent inability to fit a model based 
on the known laws of celestial mechanics to the Pioneer Doppler observations. The 
procedure can be briefly described as follows: a set of equations is used to describe 
the motion of the Pioneer spacecraft in the solar system, relative to the participating 
stations of the Deep Space Network. The motion of the spacecraft is determined by 
the initial state vector (three position and three velocity components) and perhaps 
additional parameters, such as values that characterize the velocity change as a 
result of maneuvers. Using an optimization technique such as the nonlinear least 
squares method, best fit values for the initial state vector and additional parameters 
are obtained, minimizing the difference between calculated and observed Doppler 
values. 
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When the model is modified to include an additional fictitious force that results 
in a constant sun-pointing acceleration of the sp acec raft, the statistical fit between 
the model and observations is vastly improvecPEl For this reason, it has been 
assumed that a force of unknown origin, resulting in an approximately constant, 
approximately sunward acceleration of ~ 8.74 x 10~ 10 m/s 2 , acts on the spacecraft. 

Since the initial discovery of this discrepancy between computed and observed 
Doppler readings, many attempts have been made to explain the anomaly using 
both conventional and new physics. To date, none of these attempts produced a 
satisfactory result. In particular, the possibility that the anomalous result is due 
partially or entirely to on-board forces could not be excluded; however, insufficient 
data was available to properly characterize on-board forces, resulting in wildly dif- 
fering estimates, including assertions that on-board forces are far from sufficient to 
explain the anomalous acceleration^ ^ to s uggestions that the anomalous accelera- 
tion is due entirely to on-board 

forcePM] 



3. Orbit determination 

The data that were used to develop the most precise estimate to date of the 
anomalous acceleration of Pioneer 10/1 1^ have been made available to interested 
researchers. Since the anomalous acceleration has been confirmed by several inde- 
pendent teams, it is reasonable to expect that any new attempt at orbit determi- 
nation should be able to reproduce these well established results. The anomalous 
acceleration, therefore, serves as a useful test case to validate orbit determination 
software. 

Our orbit determination is based on well established methods and techniques"^. 
It can be briefly described as follows: we integrate with high precision the (optionally 
relativistic) equations of motion for the Pioneer spacecraft in the gravitational field 
defined by major solar system bodies. The position of solar system bodies, along 
with the positions of ground-based Deep Space Network stations, are obtained from 
a software library provided by JPL. We compute the Doppler shift of the radio 
signal sent to, and received from, the spacecraft, applying all necessary corrections. 



3.1. Timekeeping 

Equations of motion are integrated using Ephemeris Time (ET) . Converting from 
ET to atomic time (TAI) can be accomplished using the following approximate 
formula (see Eqs. 2-26 through 2-28 of Ref.©: 

ET -TAI = 32.184+ 1.657 x 10" 3 sin^, (1) 

where 

E = M + 0.01671 sin M, (2) 
M = 6.239996 + 1.99096871 x 10~ 7 i, (3) 
and t is ET in seconds past J2000. 
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Data files from tracking stations are using UTC timestamps. The difference 
between UTC and TAI is an integer number of seconds ranging from 10 to 32 as 
leap seconds were added (UTC is a discontinuous time scale.) 

3.2. Equations of motion 

To calculate the orbit of the Pioneer probes, we used the Newtonian equations of 
motion without relativistic corrections. In the outer solar system, this provides a 
solution with adequate accuracy. For a spacecraft moving in the solar system that 
consists of n bodies, the spacecraft's motion is governed by the equation 

n 

where r is the spacecraft's position vector, m is the mass of the i-ih solar system 
body, multiplied by the gravitational constant, and is the position vector of the 
i-th solar system body, while the dot denotes differentiation with respect to time. 

We carry out precision integration of (|4]) using Stoermer's rule and Richardson's 
extrapolation technique. The integrator was validated by computing the orbits of 
various solar system bodies and comparing the results against data from the JPL 
HORIZONS databas<B 

3.3. Relativistic equations of motion 

In the presence of a strong gravitational field or when velocities are high, the equa- 
tions of motion (U|) must be modified, taking into account effects due to general 
relativity. Equation 4-26 in Ref. 8] takes into account all effects up to (v/c) 2 using 
the parametric post-Newtonian (PPN) formulation of general relativistic theories. 
Moyer's formulation expresses relativistic corrections in terms of the positions, ve- 
locities, and accelerations of interacting bodies. For Pioneer, we can safely omit 
terms that contain accelerations. We can also drop terms containing factors of 
1 / 1 r j — rfc| where j ^ k are indices representing solar system bodies; these terms 
represent nonlinear contributions of pairs of distant bodies, and do not significantly 
affect the Pioneer orbits. After some trivial algebra, we can express the resulting 
relativistic equations of motion in the form 

.. >^ Ui[(n - r)Ai + Bi] 



4=1 



where A4 and B 4 , which are themselves functions of the positions and velocities of 
the gravitating bodies involved and the test particle, can be written as 

2 "1 




2 +(l+ 7 )r2-2(l+ 7 )f.r 4 -^ 



(r-r,) 



(6) 



a See http://ssd.jpl.nasa.gov/7horizons 
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and 



Bi = — {(r - Ti ) ■ [(2 + 2 7 )r - (1 + 2 7 )r 4 ]} (r - U). 

n 4 



(7) 



In the non-relativistic limit (c — > oo), Ai ~ 1 and = 0, and we get back the 
Newtonian equations of motion. 

The parameters (3 and 7 are parameters of the post-Newtonian formulation; for 
Einsteinian relativity, they both have the value of 1. 



3.4. Oblateness 

Oblateness of a planetary body modifies its gravitational field. The difference be- 
tween the actual potential field of an oblate body and a spherically symmetric field 
can be described using spherical harmonics. 

The Pioneer spacecraft had close encounters with Jupiter and (for Pioneer 11 
only) Saturn. Both of these planetary bodies are rotationally symmetrical. Their 
gravitational potential field due to body i is described by the following equation: 



where Pk{x) is the fc-th Legendre polynomial in x, ai is the equatorial radius of 
planet i, 9 is the latitude of the spacecraft relative to the planet's equator, and 
is the fc-th spherical harmonic coefficient of planet i. 

In order to put ([5]) to use, first it must be translated into an expression for force 
by calculating its gradient. Second, it is also necessary to express the position of the 
spacecraft in a coordinate system that is fixed to the planet's center and equator. 

In a Cartesian coordinate system with the X1X2 plane fixed to the planet's equa- 
tor, we have 




00 





(9) 



9 = arctan 



(10) 



The gradient of the potential field can be expressed as 



dU t dUi dr dUi 89 




(11) 



dxj dr dxj 89 dxj 
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Individual partial derivatives can be calculated as 



dr Xj 
dxj r 
d9 xjx 3 



(12) 



(13) 



dxj r 2 yj{x\+x 2 2 y 



86 _ y/ X \ + xl. 
dx 3 r 2 



(14) 




(15) 



dUj _ ^ 4 a?P^(sing)cosg 

QQ ~ r k+l 



(16) 



where Pj,(x) = dPk(x)/dx. 

To use these formulae, the position of the spacecraft must be expressed in a 
coordinate system that is rotated to coincide with the planet's equatorial plane. 
The direction of the planet's North pole (which is normal to the equatorial plane) 
is known in the form of its right ascension (a) and declination (<5) angles. 

To express the coordinates of a spacecraft using a frame of reference with its 
origin at the center of planet i, one only needs to calculate r — r^. To express this 
in a frame of reference rotated by a in the X1X2 plane, we have 



Finally, to express this in a frame of reference further rotated by S in the x[x' 3 
plane, we have 



These two transformations can be combined into a single transformation matrix 
that can be used to translate coordinates to an equatorial coordinate system; the 
inverse of this matrix can be used to translate any resulting forces back into the 
original frame of reference. 

3.5. Time delay 

Our observable is radio-metric data, measured at specific moments in time at specific 
Earth stations. To accurately compute a matching orbit, we need to know not where 
the spacecraft was at the moment the signal was received on the Earth, but where 
the spacecraft was at the moment it transmitted that signal. 



x\ = x\ cos a + X2 sin a, 
x' 2 = —xi sin a + x 2 cos a, 
x' 3 = x 3 . 



(17) 
(18) 
(19) 




x' x sin 8 — x 3 cos S, 
x\ cos (5 + ^3 sin S. 



(20) 
(21) 
(22) 
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The time delay can be calculated using an iterative process that is known to 
converge rapidly. We assume that the signal propagates from the transmitter at 
Pi(t) to the receiver at P2(£)- If the time of reception ti is known, the location 
of the receiver can be accurately calculated as ~P%{t'i). We use ti as a first crude 
estimate for the time of transmission: 

*i(o) = h- (23) 

The corresponding location along the transmitter's path is Pi(t 1 ( )), which we use 
to estimate the propagation delay and arrive at an improved estimate of the time 
of transmission: tux) = t% — |Pi(^i(o)) ~ ^2{t%)\/ c. This process can be improved 
iteratively, using 

*i« = *2 - ~|Pi(*i(i-i)) - P 2 (* 2 )|. (24) 

By itself, this procedure alone does not provide a sufficiently accurate estimate 
for the time delay. The reason is that the radio signal of the Pioneer probes travels 
through the gravitational field of the Sun and the planets, where the effective speed 
of light is less than c. The consequent delay At, known as the Shapiro time delay, 
between two points Pi and P2 can be calculated as (cf. Ref. [8]Eq. 8-54): 

At = lln ri+r2+ri2+ ; , (25) 

where / = //(l + 7)/c 2 ; /1 is the mass of the body whose potential field is being 
considered, multiplied by the gravitational constant; 7 is a parameter of the PPN 
(Parameterized Post-Newtonian) formalism that is equal to 1 for general relativity; 
r\ and T2 are the distances to Pi and P2 from the gravitating body; and r*i2 = 
Pi — P2I is the distance between Pi and P2. 

For a spacecraft moving in deep space far from planetary bodies, only the Sun's 
gravitational potential needs to be considered for the Shapiro delay. Furthermore, 
the distance of the spacecraft from the solar system barycenter can be substituted 
in place of the distance from the Sun's center of gravity, resulting in a substantial 
simplification of the calculations at no cost in terms of accuracy. 

3.6. Calculating the Doppler observable 

Our observable is Doppler: that is to say, the goal of the orbit determination exercise 
is to find an initial state vector that minimizes the difference between the measured 
frequency of the spacecraft's radio transmissions vs. the calculated frequency. As the 
spacecraft lacked a precision oscillator on board, this always means two-way or three- 
way Doppler: a signal of known frequency is transmitted to the spacecraft, which 
in turn reradiates a coherent response signal to the Earth, that is then received 
by either the same antenna that was used to do the transmission, or a different 
antenna. 

If we know the exact positions and velocities of the transmitting and receiving 
stations as well as the spacecraft itself, the frequency of the received signal can 
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be calculated to high precision. In an inertial coordinate system, with the effects 
of special relativity taken into account, a signal of frequency f\ transmitted from 
station Pi (moving with velocity vi) will be received at station P2 (moving with 
velocity V2) at the frequency 



Doppler is not an instantaneous measurement: the receiving station counts the 
difference in cycles between the received frequency and the reference frequency for 
a specific interval. When the spacecraft is in deep space and the forces acting on 
it are very small, we can use the spacecraft's average velocity during the Doppler 
interval in ([25)) . However, this can lead to significant errors when the spacecraft is 
subjected to significant acceleration, e.g., when it is near a major solar system body. 

Note that for the Earth stations, their instantaneous velocity must be accounted 
for properly, as the Earth stations are always accelerating as a result of the Earth's 
rotation and orbital motion. 

3.7. Nongravitational forces and effects 

In addition to gravity, several small forces act on the spacecraft and/or the trans- 
mitted radio signal that must be taken into account. 

3.7.1. Solar pressure 

Light from the Sun produces a noticeable amount of force on the spacecraft. Near 
the Earth, every square meter of surface receives 1366 W of solar power ("solar 
constant"); the total area of the Pioneer HGA is approximately 5.91 m 2 , translating 
into more than 8 kW! (In contrast, the anomalous acceleration of the Pioneers is 
equivalent to ~65 W of power emitted in a collimated beam in the appropriate 
direction.) Far from the Sun, solar pressure is significantly lower, but even at 20 AU, 
the HGA receives ~20 W of sunlight, and therefore this effect cannot be ignored. 

Since most of the time, the HGA is pointed towards the Earth and, consequently, 
in the approximate direction of the Sun, we need not model the spacecraft body 
behind the HGA; instead, we can use a "flat disk" model of the HGA. The HGA 
surface area is known. Sunlight that reaches the surface will transfer its momentum 
to the spacecraft; in addition, sunlight that is reflected off the surface, or reemitted 
by the surface, transfers additional momentum. 

It is obvious that three principal directions need to be considered. Intercepted 
light will transfer momentum that is in the direction of the Sun-spacecraft line. 
Diffusely reflected, as well as re-emitted solar energy will produce momentum along 



(26) 




(27) 
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the spin axis. Lastly, specularly reflected sunlight will produce momentum along 
a line that is the reflection of the Sun-spacecraft line along the spin axis, in the 
Sun-spacecraft-spin axis plane. 

Quantitatively, force acting on the spacecraft due to intercepted sunlight can be 
calculated as 

„ sA(r ■ n) , . 

Fintcpt = V 4 ' r, 28 

where s is the normalized solar constant (s = 1366 W/m 2 x(l AU) 2 ), A = 5.91 m 2 
is the surface area of the HGA, n is a normal vector in the direction of the spin 
axis pointing towards the rear of the spacecraft, and r is the Sun-spacecraft vector. 
As the value of A is somewhat uncertain (the antenna is not exactly round, it has 
cutouts, and parts of it are shadowed by structural elements and the medium gain 
antenna) we may multiply this expression with a proportionality factor a ~ 1: 

„ sA(r ■ n) 

Fintcpt = a y j ; r, 29 

The amount of specularly reflected sunlight will be proportional to the amount of 
intercepted sunlight, with a proportionality factor a. Its direction can be computed 
using elementary vector geometry: 

F sp cc = — v [Fintcpt — 2(F; n tcpt ■ n)n] . (30) 

Whatever solar radiation was not specularly reflected is either reflected diffusely, 
or it is absorbed by the HGA and reradiated. Either way, the resulting force will be 
in the direction of the spin axis. The force will be proportional to e(-Fi n tcpt — -Fspcc), 
where e is a proportionality factor that is a function of the antenna's absorptance 
and front and rear infrared emittance: 

2 
3 

Using values from Pioneer project documentation^ESIl]] we can calculate 



diff — o e (|Fintcptl — |F spoc |)n. (31) 



£ = (1 - a S r) + Osrfcfront - e rea r)/ (cfront + e rcar) ^0.98, 

where a sr ~ 0.21 is the solar absorptance of the antenna front side, ef ron t — 0.85 is 
its infrared emittance, while e rca r — 0.04 is the infrared emittance of the antenna 
rear side. The additional factor of 2/3 is characteristic of the force due to diffuse 
omittance by a flat emitter. 

The parameters a ^ 1, — 1 < e ~ 0.98 < 1, and < a < 1 can be solved 
for as an orbital solution is obtained. Alte rnativ ely, their values can be calculated 
approximately from project documentatio n^ 1 H 



3.7.2. Media effects 

The atmosphere and solar plasma can introduce a significant additional signal prop- 
agation delay. One way to deal with this delay is to ignore all data points measured 
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when the spacecraft was less than 5 degrees above the horizon as seen from either 
the transmitting or the receiving antenna, and when the angular separation between 
the Sun and the spacecraft was less than 5 degrees. However, it is also possible to 
estimate these delays. 

• Solar Plasma 

Delay due to solar plasma is a function of the electron density in the plasma. 
Although this can vary significantly as a result of solar activity, the propaga- 

I Vol 

tion delay AT (in microseconds) can be approximated using the formula : 

AT = 1.3446 x 10~ 19 x /~ 2 x J N e dl, (32) 

where / is the signal frequency, and N e is the electron density, which is 
integrated along the propagation path I. For N e , an empirical formula is 
available^ 

N e = 2.21 x 10 14 (j^j + 1.55 x 10 12 , (33) 

where Rq = 6.96 x 10 8 m is the solar radius, and r is the distance from 
the Sun along the propagation path. Although (1B"2"1) cannot be expressed in 
closed form, numerical integration can yield a table of values that can be 
used, along with an interpolation subroutine, to obtain values of AT, as a 
function of the spacecraft-to-Earth distance and the Sun-Earth-spacecraft 
angle, with sufficient precision. 

• Tropospheric Delay 

ChacU^ estimates the delay due to signal propagation through the tropo- 
sphere using the following formula: 

R = . = 1 A , (34) 

where R is the additional propagation path, E is the elevation angle, and 
A = Adry + A wct and B — B^y + B wct are parameters defined as 



Adry 


= 0.00143, 


(35) 


A W et 


= 0.00035, 


(36) 


-Bdry 


= 0.0445, 


(37) 


B wc t 


= 0.017. 


(38) 



3.7.3. Spacecraft spin 

The radio signal used to communicate with the Pioneer spacecraft is circularly 
polarized. As the spacecraft itself is spinning, every rotation of the spacecraft effec- 
tively adds an extra period to the (received or transmitted) signal. In other words, 
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the spacecraft's spin frequency (l/60th its rate of revolution in units of rpm) must 
be added to both the received and the transmitted signal. 

The sign of the spin correction can be unambiguously determined by comparing 
archived Doppler data files, some of which were corrected for spin and some of 
which weren't. This comparison clearly reveals that the spin of the spacecraft must 
be added to the frequency of both the uplink and the downlink signal with a positive 
sign. 

The spacecraft spin can be directly obtained from telemetr y^ -^ l A star sensor, 
or one of two sun sensors were used to measure the spin period of the spacecraft 
on board, and this value was telemetered to the Earth using a clock resolution 
of 1/8192 seconds. This information was sent to the Earth in the form of three 
telemetry words in the telemetry data streairM 

The sun sensor on the spacecraft was operational only up to a distance of 35 
AU from the Sun. On Pioneer 10, the star sensor failed after Jupiter encounter. As 
a result, once the spacecraft passed beyond 35 AU, the spin period could no longer 
be determined on board. Readings from the Imaging Photo-Polarimeter instrument 
were used to occasionally derive further estimates of the spin rate. 



3.7.4. Maneuvers 

After planetary encounters, the Pioneer spacecraft no longer executed trajectory 
correction maneuvers. However, they still periodically performed precession maneu- 
vers, designed to ensure that the spacecraft remains pointed in the direction of the 
Earth with sufficient accuracy to maintain reliable two-way radio communication. 

These maneuvers involved the firing of thrusters on opposite sides of the space- 
craft, one in the fore, the other in the aft direction, for equal duration. In principle, 
such a maneuver imparts no net change in momentum on the spacecraft. In prac- 
tice, uncertainties in thruster alignment and firing duration result in a small net 
change in the spacecraft velocity in a random direction. As the spacecraft spins, 
random changes in momentum in a direction perpendicular to the spin axis cancel 
out over time, only increasing the uncertainty in the spacecraft's lateral position 
in the sky slightly. However, the firings may increase or decrease the spacecraft's 
velocity along its spin axis, and this must be taken into account. 

The magnitude of this velocity change is essentially random; it cannot be pre- 
dicted. However, the exact time of precession maneuvers is known from the telemetry 
data stream. 

Specifically, the telemetry data stream provides thruster pulse countJ^^ 6-bit 
counters for each thruster that are incremented every time the thruster in question 
is fired. 

We were able to compile this information from the telemetry into a list of times 
when precession maneuvers occurred. This list is directly imported into our orbit 
determination code if an estimate of maneuvers is requested. 
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3.7.5. Heat 

The study of on-board generated heat and its effe cts on th e orbit of Pioneer 10 and 



code has the ability to deal with heat sources and calculate their contribution to 
acceleration. This capability will be used with the newly compiled Doppler data set; 
results will be reported as they become available. 

3.7.6. Radio beam recoil force 

The spacecraft's radio beam also carries momentum that corresponds to the ap- 
proximately 8 W of transmitted power. The spacecraft's HGA is always aimed at 
the approximate direction of the Earth, deviating from that orientation by no more 
than a couple of degrees. However, it has been argued-*- 1 that not all of the radio 
power emitted by the transmitter hits the reflecting dish of the HGA; some of it, 
typically 10%, would miss the antenna and would be radiated in the opposite direc- 
tion, at an approximately 45° angle. This means that the radio beam will translate 
into momentum by an approximate efficiency of ~ 83%^. 

The actual power of the transmitter was always assumed to be constant, but 
telemetry suggest s othe rwise: the transmitter power is directly measured, and it 
varied significantl y 1 '. 

3.7.7. Solar system data 

Rather than integrating the equations of motion for all solar system bodies, we 
use the planetary ephemerides published by JPL. Specifically, we are using JPL's 
CSPICE libranM to obtain high accuracy planetary positions. To calculate the orbit 
of the Pioneer spacecraft, we utilize the state vectors of the Sun and the eight major 
planets from the DE-414 ephemerides. 

Except for cases when the spacecraft are near Jupiter or Saturn, we ignore the 
effects of individual moons. We also ignore smaller bodies such as Pluto, Ceres, 
Charon, or the larger Kuiper Belt objects. 

The CSPICE library is also used to obtain accurate estimates of planetary 
masses, multiplied by the gravitational constant. Our work therefore does not de- 
pend on the accurate value of the gravitational constant, and all calculations are 
based on data sets with a consistent origin. 

We also used the CSPICE library to obtain accurate state vectors of the par- 
ticipating stations of the Deep Space Network. The library allows us to obtain 
momentary positions and velocities. However, a typical Doppler measurement has a 
duration ranging from several seconds to several minutes, during which the station 
moves significantly in a solar system barycentric frame of reference, and its velocity 
changes as well. In a geocentric reference frame, the station moves along a circular 



11 are a major focus of our on-going 




ion 



b See http : //naif . jpl . nasa . gov/ 
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arc in the plane of the equator. For an arc spanning an angle <f> with radius r, the 
center of gravity is located at angle <p/2 at a distance of [(sin0/2)/(</>/2)]r (</> mea- 
sured in radians) from the origin. This expression was used to correctly average the 
station position due to the Earth's rotation. 

The right ascension and declination of Jupiter's and Saturn's north pole are pro- 
vided in the appropriate CSPICE kernels, along with spherical harmonic coefficients 
(J 2 through J 8 ). 

In the vicinity of the Earth, Jupiter, and Saturn, the gravitational influences of 
individual moons must be considered. Systems of satellites need not be considered 
when the spacecraft is far from the planetary system. A distance of 10 11 m (~ 2/3 
AU) works as a sensible cut-off. 

Jovian and Saturnian satellite ephemerides are obtained from the CSPICE kernel 
files JUP204 (Galilean satellites, DE-405), JUP256 (minor satellites, DE-405), and 
SAT252 (Saturnian satellites, DE-414.) 

Station position and velocity information (in the Earth-centered ITRF93 refer- 
ence frame) for all operating DSN stations are readily available in CSPICE kernel 
files. It is somewhat harder to obtain station positions for DSN stations that have 
been decommissioned, or DSN stations that were relocated. 

Specifically, position and velocity information for DSS-11, DSS-41, DSS-44, DSS- 
51, and DSS-62 were provided by FolkneJl, along with position information for 
DSS-12 and DSS-61 before their relocation in 1978-79. An ambiguity concerning 
the exact date of the relocation of DSS-12 was resolved by Ref. [TrJl 

3.7.8. Orbit data files 

Actual measurements of the Pioneer spacecraft's radio frequency are stored in files 
called ODFs (Orbit Data Files.) These files have a simple binary format encoding 
the following pieces of information: 

• Time of reception 

• Count period 

• Reference frequency 

• Transmitting station identifier 

• Receiving station identifier 

• Doppler observable. 

In principle, the Doppler observable can be reconstructed from the first five of 
these values: from the time of reception and the known orbit of the spacecraft, taking 
light trip time into account, the time of transmission on board the spacecraft can 
be calculated, and from that, the time of transmission at the transmitting station 
can be obtained. The relative velocities of these three locations, combined with the 



c Private communication. 
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reference frequency, can be used to obtain the received frequency; the difference 
between this and the reference frequency is the Doppler reading. 

One complicating factor is that the actual frequency of transmission might not 
have been the reference frequency; indeed, the reference frequency itself might have 
been changing in time as well, i.e., in many cases these frequencies were "ramped". 
The ODFs contain ramp data records that can be used to determine the precise 
transmission and reference frequencies at specific DSN stations at any given time. 
The absence of a ramp record typically manifests itself as an offset of several kHz 
between the calculated and observed frequency; Doppler data records with missing 
ramp information can thus be easily identified and rejected. 

The ODFs also had other correctable issues, such as incorrect encoding of fre- 
quencies (e.g., marking an X-band frequency as a Ku-band frequency) that were 
easily fixed. 

A far more insidious problem is related to spin correction. Often, ODFs are ad- 
justed to account for the spacacraft's spin, but this is not always the case. Notably, 
the 11.5 year Pioneer data set contained records both with and without spin correc- 
tion; data points prior to July 18, 1990 had spin correction applied, whereas later 
data points were uncorrected^. 



4. Results 

The considerations presented so far were used to build an orbit determination pro- 
gram that took into account solar system ephemerides published by JPL and uti- 
lized through the CSPICE library, DSN station locations on the Earth, again cal- 
culated using the CSPICE library and averaged for the duration of Doppler inter- 
vals, Doppler observations obtained from ODFs, spacecraft spin data obtained from 
telemetry, the spacecraft's trajectory, light trip times including the Shapiro delay, 
a flat disk solar pressure model, and optionally, additional telemetry /design values 
used to estimate thrust due to heat and the radio beam recoil force. The goal of 
this application was to minimize the sum of the squares of differences between com- 
puted and observed Doppler; the minimization algorithm used was the nonlinear 
least squares algorithm with approximate Jacobian from the standard MINPACK 
librarjo We weighted individual data points using the Doppler duration, with the 
assumption that the standard deviation of the data is the inverse of that duration; 
for instance, if the Doppler count was 1000 seconds, we assume that the data point 
is accurate to 1 ms within 1 standard deviation. 

With all the tools in place, the first order of business was to use these to confirm 
the earlier results regarding the Pioneer anomaly. 



d Private communication from P. Tortora. 
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4.1. Confirmation of the Pioneer anomaly 

To confirm previous results regarding the anomalous acceleration of Pioneers 10 and 
11, we used ODFs provided by JPL. These ODFs contained verified and filtered 
records spanning the time period 1986-1998 (Pioneer 10) and 1986-1990 (Pioneer 
11). These were the exact same data records examined by JPL in 

For the initial state vector, we used the approximate Pioneer ephemerides from 
JPL HORIZONS. This produced a good fit already, with the Doppler residual for 
most data points not exceeding a few Hz (Figure [2] top left and Figure [3] top left.) 




Fig. 2. Pioneer 10 Doppler residuals. Vertical axis is residual frequency in Hz; horizontal axis is 
time (in years). Dotted vertical lines correspond to maneuvers. Top left: residuals prior to fitting. 
Top right: residuals after a best fit of the initial state vector, the anomalous acceleration and 28 
maneuvers. Bottom left: residuals after the anomalous acceleration has been removed from the 
best fit shown in top left. Bottom right: residuals after a best fit of the initial state vector and 28 
maneuvers, without anomalous acceleration. 



4.2. Pioneer 10 results 

For Pioneer 10, we used JD (Julian Day) 2446765 through JD 2451016 as the period 
of investigation; the data set contained 19758 Doppler data points in total. About 
one third of the Doppler data points were obvious outliers (yielding residuals of 
up to several MHz). Most of these data points were in the last two years of the 
mission, corresponding to changes in DSN station configuration, and therefore, it 
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Fig. 3. Pioneer 11 Doppler residuals. Vertical axis is residual frequency in Hz; horizontal axis is 
time (in years). Dotted vertical lines correspond to maneuvers. Top left: residuals prior to fitting. 
Top right: residuals after a best fit of the initial state vector, the anomalous acceleration and 30 
maneuvers. Bottom left: residuals after the anomalous acceleration has been removed from the 
best fit shown in top left. Bottom right: residuals after a best fit of the initial state vector and 30 
maneuvers, without anomalous acceleration. 



is possible that these ODF records are misinterpreted by our code. Nevertheless, 
after the removal of outliers, a sufficient number of good data points remains for 
analysis. The time period included 28 maneuvers. We used an incremental filtering 
strategy to process this data set. We used a nominal spacecraft mass of 250 kg, 
corresponding to a loss of one quarter of the propellant on board, when calculating 
nongravitational forces, such as solar pressure. In addition to filtering outliers, we 
also set a minimum horizon angle of 5° and a minimum angular separation angle of 
5° between the spacecraft and the Sun. Our initial run, prior to parameter fitting, 
yielded 13525 data points. 

To fit the data, we followed an incremental strategy. First, we estimated the 
initial state vector, maneuvers, and the constant sunward acceleration using every 
20th record in the data set. This allowed for a rapid process to find good initial 
estimates for these parameters, allowing us to iteratively reduce the filter bandwidth; 
eventually, we were able to narrow the bandwith to 100 mHz, rejecting data points 
that yielded a Doppler residual in excess of 50 mHz. At the end of this stage, we 
were left with 13119 Doppler data points. 

Finally, we fitted these data points without any additional filtering or sampling. 
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We obtained a result of (9.03±0.13) x 10 -10 m/s 2 . Combined with error terms from 
Ref. [5] (but excluding terms related to thermal modeling, which is the subject of 
our on-going effort), amounting to ±0.85 x 10 -10 m/s 2 , we obtain 



with a root mean square (rms) Doppler residual of 10.1 mHz (Figure [2] top right). 

Figure [5] (bottom left) shows the Doppler residuals using the fitted value of 
the initial state vector and maneuvers, but with ap set to 0. The result is an 
unambiguous drift, another clear indication of the anomalous acceleration. 

Is it possible to fit the data with ap = 0? Clearly, much of the anomalous 
behavior can be absorbed into the 28 maneuvers that took place during this time 
period. However, one would expect a worse fit in this case. Indeed, this is what we 
found; after fitting 13111 data points by varying only the initial state vector and 
maneuvers, but without an anomalous acceleration component, the rms residual 
was 43.8 mHz (Figure [5] bottom right). This clearly indicates that the constant 
acceleration term offers a statistically significant, superior fit to the Doppler data. 

4.3. Pioneer 11 results 

The results for Pioneer 11 are similar to the Pioneer 10 results. The Doppler data, 
from JD 2446431 to JD 2448256, covered approximately five years, with a total 
of 11423 data points, containing only a small number of outliers. The nominal 
spacecraft mass used was 232 kg, to account for the increased use of propellant 
(three quarters of the propellant on board) during Pioneer ll's trajectory correction 
maneuvers. For Pioneer 11, a shorter data span was available, but it contained more 
data points; there were also more maneuvers performed during the study period. 
Our initial filtering yielded 11230 data points (Figure [3] top left). 

Following an incremental editing strategy just like in the case of Pioneer 10, we 
were able to fit the initial state vector, 30 maneuvers, and the anomalous accelera- 
tion. With 11080 data points remaining in the data set after filtering, we obtained 
(8.21 ± 0.66) x 10~ 10 m/s 2 , which, combined with error terms from Ref. 2, yields 



with an rms residual of 8.2 mHz. The substantially larger formal error can be 
attributed to the shorter timespan covered by the data set, and more frequent 
maneuvers (Figure [3] top right). 

Figure [3] (bottom left) shows the computed Doppler residual using the fitted 
initial state vector and maneuvers, but with ap set to zero. Again, the anomalous 
velocity drift is clearly indicated in the data. 

As in the case of Pioneer 10, we also attempted to fit the Doppler data using 
the initial state vector and maneuvers only, with the anomalous acceleration set to 
zero. The result, shown in Figure [3] (bottom right), is inferior, though not as clearly 
so as was the case for Pioneer 10. The rms residual remains a modest 12.4 mHz on 



apio = (9.03 ± 0.86) x 10~ 10 m/s 2 , 



(39) 



apn = (8.21 ± 1.07) x 10~ 10 m/s 2 , 



(40) 
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11070 data points, due to the fact that the shorter timespan and greater frequency 
of maneuvers allowed a tighter fit even without ap. 

4.4. Spectral Analysis 

How good are the results that we obtained so far? As can be seen from the unfitted 
residuals (Figure top left and Figure [3] top left), any mismodeling is bound to 
produce a clear annual signature, due to the orbital motion and rotation of the 
Earth where the DSN antennae are located. 




Fig. 4. Pioneer 10 residuals normalized spectral density. Left: residuals prior to fitting; a strong 
annual signal (position indicated by vertical green line) is clearly seen. Right: residuals after the 
fit depicted in Figure [2] (top right). The annual signature has vanished, and only a small peak is 
visible at a period of 1 day, indicating a small diurnal mismodeling. 




Fig. 5. Pioneer 11 residuals normalized spectral density. Left: residuals prior to fitting; a strong 
peak is visible at a period of ~370 days, which is an annual signal the period of which is skewed 
by the spacecraft's proper motion. Right: residuals after the fit depicted in Figure \3\ (top right). 
No obvious peaks are visible that would indicate a diurnal or annual mismodeling. 



Such signatures must be observable in the power spectrum of Doppler residuals. 
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Because the sampling intervals are not uniform (there are dropouts and blackouts 
in the Doppler data stream) , traditional FFT methods are not readily applicable to 
our resi duals; however, the residuals can be readily analyzed using Lomb-Scargle 
methodsP^lMI 

The results for the unfitted data sets, shown in Figure[4](left) and Figure [5] (left) 
clearly indicate the presence of an annual signal (with a period slightly longer than 
one year due to the spacecrafts' lateral motion). A small peak is also visible at a 
period of 1 day, though it is not significant compared to other peaks. 

After the fit, the annual terms have all but vanished from the power spectrum 
(Figure [4] right, and Figure [5] right) . For Pioneer 10, a noticeable peak is present at 
a period of 1 day, indicating a possible diurnal mismodeling; no such peak can be 
observed in the Pioneer 11 post-fit power spectrum. 

4.5. Jerk term 

A recurring question concerning the Pioneer anomaly is its constancy. To test the 
extent to which the anomalous acceleration is constant in time, we implemented the 
ability to estimate a secondary acceleration, i.e., "jerk" term in our orbital solution. 
For Pioneer 10, we obtained 

api = (10.96 ± 0.89) x 10~ 10 m/s 2 , (41) 
a P10 = (-0.21 ± 0.04) x 10~ 10 m/s 2 /year, (42) 

with an rms residual of 9.8 mHz on 13130 data points. 
For Pioneer 11, the result was 

opu = (9.40 ± 1.12) x 10~ 10 m/s 2 , (43) 
dpn = (-0.34 ± 0.12) x 10~ 10 m/s 2 /year, (44) 

with an rms residual of 8.1 mHz on 10995 data points. 

We note that this jerk term is consistent with the expected temporal variation 
of a recoil force due to heat generated on board and emitted anisotropically. 

5. Plans and conclusions 

The results presented in this work confirm the existence of the Pioneer anomaly. 
Doppler observations from the twin Pioneer spacecraft only match calculated orbits 
if we assume the presence of a small, approximately constant anomalous acceleration 
in the sunward direction. 

We have demonstrated our ability to reproduce results using independently de- 
veloped orbit determination software with sufficient accuracy. In addition to the 
Doppler data files, our software has the ability to process telemetry data obtained 
from the spacecraft. Specifically, telemetry records are used to establish the space- 
craft's spin and the timing of maneuvers. The telemetry records can also be used to 
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estimate thermal emissions from the spacecraft, which is the subject of our on-going 
study. 

In addition to confirming the numerical value of the anomalous acceleration, 
we also considered a secondary acceleration ("jerk") term. We demonstrated that a 
moderate jerk term is consistent with the Doppler data, and therefore, an anomalous 
acceleration that is a slowly changing function of time cannot be excluded at present. 

Our modeling can be improved. We can presently model planetary flybys with 
limited precision; the maximum Doppler residual during flybys of Jupiter and Sat- 
urn is as high as a few Hz (indicating a velocity mismodeling of several 10 cm/s) 
leaving room for improvement. We also plan to make our tropospheric model more 
sophisticated; even in the absence of historical weather records, Markwardt sug- 
gests that seasonal weather data can be usedl to achieve good modeling accuracy. 
We may also model signal propagation delays due to charged particles in the Earth's 
ionosphere. 

Our main effort, however, is aimed at developing robust orbit estimates that 
incorporate a telemetry-driven time dependent model of on-board generated thermal 
recoil forces. The goal is to determine unambiguously the extent to which thermal 
forces may be responsible for the anomalous acceleration of the Pioneer spacecraft. 

Acknowledgments 

I would like to thank Dr. Slava Turyshev of JPL for his support, advice, and en- 
couragement; and Larry Kellogg, formerly of NASA Ames (now retired), who sin- 
glehandedly saved the telemetry records of Pioneer 10 and 11. 

References 

1. John D. Anderson, Philip A. Laing, Eunice L. Lau, Anthony S. Liu, Michael Martin 
Nieto, and Slava G. Turyshev. Indication, from Pioneer 10/11, Galileo, and Ulysses 
Data, of an Apparent Anomalous, Weak, Long- Range Acceleration. Phys. Rev. Lett., 
81:2858-2861, 1998. 

2. John D. Anderson, Philip A. Laing, Eunice L. Lau, Anthony S. Liu, Michael Martin 
Nieto, and Slava G. Turyshev. Study of the anomalous acceleration of Pioneer 10 and 
11. Phys. Rev. D, 65:082004/1-50, 2002. 

3. Slava G. Turyshev, Viktor T. Toth, Larry R. Kellogg, Eunice L. Lau, and Kyong J. Lee. 
The Study of the Pioneer Anomaly: New Data and Objectives for New Investigation. 
Int. J. Mod. Phys. D, 15:1-55, 2006. 

4. Craig B. Markwardt. Independent Confirmation of the Pioneer 10 Anomalous Accel- 
eration. 2005. 

5. J. I. Katz. Comment on "Indication, from Pioneer 10/11, Galileo, and Ulysses Data, 
of an Apparent Anomalous, Weak, Long-Range Acceleration". 1998. 

6. Edward M. Murphy. A prosaic explanation for the anomalous accelerations seen in 
distant spacecraft. 1998. 

7. Louis K. Scheffer. Conventional Forces can Explain the Anomalous Acceleration of 
Pioneer 10. Phys. Rev. D, 67, 2003. 



Private communication. 



23, 2009 1:2 WSPC/INSTRUCTION FILE pioneer orbit 



22 

8. Theodore D. Moyer. Formulation for Observed and Computed Values of Deep Space 
Network Data Types for Navigation. 2000. 

9. Oliver Montenbruck and Bernhard Gill. Satellite Orbits. Springer, third edition, 2005. 

10. NASA Ames Research Center. Pioneer F/G: Spacecraft Operational Characteristics. 
1971. PC-202. 

11. TRW Systems Group. Pioneer F/G: Thermal Control Subsystem Design Review Num- 
ber 3. 1971. 

12. Ho, C. and Morabito, D. Solar corona and solar wind effects. In DSMS Telecommuni- 
cations Link Design Handbook 106, Rev. A., 2005. 

13. Sovers, O. J. and Jacobs, C. S. Observation Model and Parameter Partials for the JPL 
VLBI Parameter Estimation Software "MODEST" -1994. JPL Publication, (83-39), 
1994. 

14. Viktor T. Toth and Slava G. Turyshev. The Pioneer Anomaly: Seeking an explanation 
in newly recovered data. Can. J. Phys., 84:1063-1087, 2007. 

15. V. T. Toth and S. G. Turyshev. Pioneer Anomaly: Evaluating Newly Recovered Data. 
In AIP Conf. Proc, volume 977, pages 264-283, 2008. 

16. Douglas J. Mudgway. Uplink-Downlink: A History of the Deep Space Network 1957- 
1991. 2001. NASA SP-2001-4227. 

17. William H Press, Saul A. Teukolsky, William T. Vetterling, and Brian P. Flannery. 
Numerical Recipes in C, Second Edition. Cambridge University Press, 1992. 

18. Earl F. Glynn, Jie Chen, and Arcady R. Mushegian. Detecting periodic patterns 
in unevenly spaced gene expression time series using Lomb-Scargle periodograms. 
Bioinformatics, 22(3):310-316, 2006. 



